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Conditional Eulerian and Lagrangian velocity increment statistics of fully 
developed turbulent flow 
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Conditional statistics of homogeneous isotropic turbulent flow is investigated by means of high-Reynolds 
number direct numerical simulations performed with 2048 3 collocation points. Eulerian as well as Lagrangian 
velocity increment statistics under several conditions are analyzed and compared. In agreement with ex- 
perimental data longitudinal probability density functions P(S^u\ei) conditioned on a scale-averaged energy 
dissipation rate are close to Gaussian distributions over all scales within the inertial range of scales. Also trans- 
verse increments conditioned on either the dissipation rate or the square of the vorticity have quasi-Gaussian 
probability distribution functions (PDFs). Concerning Lagrangian statistics we found that conditioning on 
a trajectory averaged energy-dissipation rate e T significantly reduces the scale dependence of the increment 
PDFs P(S T Ui\e T ). By means of dimensional arguments we propose a novel condition for Lagrangian incre- 
ments which is shown to reduce even more the flatness of the corresponding PDFs and thus intermittency in 
the inertial range of scales. The conditioned Lagrangian PDF corresponding to the smallest increment con- 
sidered is reasonably well described by the K41-prediction of the PDF of acceleration. Conditioned structure 
functions show approximately K41-scaling with a larger scaling range than the unconditioned ones. 

Keywords: Homogeneous isotropic turbulence, conditional statistics, intermittency 



I. INTRODUCTION 

The problem of anomalous scaling can be seen as 
one of the great unsolved problems in turbulence re- 
search. The scaling laws of velocity structure functions 
S p (l) — ((5iu) p ) ~ fa within the inertial range of scales 
of fully developed turbulence has inspired a lot of pub- 
lications over the last decades^ 2 -. The velocity incre- 
ment under consideration is usually either the longitu- 
dinal Sfu = (u(x + I) — u{x)) ■ 1 or the transverse one 
5\ ; u = \(u(x + I) — u{x)) x l\. Both are so called Eule- 
rian increments because the velocity differences are taken 
over spatial separations at a the same instant of time. 
Kolmogorov's K41-theory2 implies a linear scaling law 
Cp = p/3 not distinguishing between the two different 
types of increments mentioned before. However, direct 
numerical simulations (DNS) and experiments show a de- 
viation of the formi 



of election-MHD turbulence^ the differences between lon- 
gitudinal and transverse scaling exponents were found to 
decrease with Reynolds-number. One has also to be very 
careful in the determination of these scaling exponents 8 
as longitudinal and transverse structure functions possess 
differing scaling ranges. 

Guided by Kolmogorov's refined self similarity 
hypothesis^ (RSH) 



(2) 



which states a relation between the local energy dissipa- 
tion rate 



e(x) 



■diUj(x)f 
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averaged over a scale I 



(1) 



e(x + sl)ds 



(4) 



with a positive /i p . The question whether longitudinal 
and transverse statistics possesses two different sets of 
scaling exponents (p , Cp respectively, is still under dis- 
cussion. Experimental observations^ as well as DNS 6 
found slightly smaller scaling exponents for the high- 
order transverse than for the longitudinal structure func- 
tions. It is not yet clear whether these findings are finite- 
Reynolds number and/or anisotropy effects. In the case 



(I indicating the same line appearing in s)u and $i~u) 

and the velocity fluctuation S\ u over that scale Gagne et 
Ac— experimentally measured conditional velocity in- 
crement statistics. They found the probability density 
functions (PDFs) P(5^u\ei) to be nearly Gaussian from 
the dissipation- up to the integral-scale which implies a 
linear scaling law with \i v = in ([lj. It is believed that 
anomalous scaling (/x p ^ 0) in turbulent flows has its ori- 
gin in small-scale intermittency of the local energy dis- 
sipation rate. This point of view is supported by this 
result, namely that the statistics of increments become 
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Gaussian once they are conditioned on a scale-averaged 
energy dissipation rate. 

Recently, new experimental techniques have provoked 
a renewed interest in Lagrangian statistics^— . Here, 
velocity increments are taken along trajectories of fluid 
elements (tracers). Lagrangian velocity increments are 
defined by 

S T Vi = Vi(r) - Vi(0) 

= Ui(X(x o ,T),T)- Ui (X{x ,0),0), [ ' 

where X(xq, t) denotes the trajectory of a tracer which 
started at the position Xq at time t — 0. Although one 
might expect a scaling law of the corresponding struc- 
ture functions S^{1) — ((S T Vi) p ) ~ 7-C* within the tem- 
poral inertial range of scales it has not yet been clearly 
observe d 14 ! 15 . The origin of the strong measured inter- 
mittency, the very existence of an inertial range of scales, 
the corresponding scaling exponents, and their relation 
to Eulerian intermittency are still open issues 1 ^—. 

In this paper we measure conditional velocity statistics 
both in the Eulerian as well as in the Lagrangian frame 
of reference by means of high-Reynolds number DNS. We 
approve the experimental results obtained by Gagne et 
alJ£ and Naert et alJ£ and complement their findings 
by a detailed scale by scale analysis and an investigation 
of the statistics of conditioned transverse velocity incre- 
ments. Furthermore we analyze conditioned Lagrangian 
increment statistics. 

In the next section we briefly present the numerical 
method. Section Unl presents the results in the Eulerian 
frame and section IIVI those in the Lagrangian frame of 
reference. Conclusion are summarized in section [Vj 

II. NUMERICS 

The numerical simulations were performed by solving 
the incompressible Navier-Stokes equations 

d t u + (u ■ V)« = / - Vp + vAu (6) 
V • u = 0, (7) 

in a periodic cube with a pseudo-spectral method using 
a high-order exponential cut-of f 20 ! 21 . 

We parallelize the computations via a pencil geome- 
try by means of the San Diego P3D-FFT— and explore 
the BlueGene/P-architecture (the 2048 3 simulation was 
performed on 32k processors). The time integration of 
the velocity field is done by means of a strongly stable 
Runge-Kutta third order scheme^. In order to maintain 
a statistically stationary flow a forcing / is applied which 
keeps constant the modes of the two lowest Fourier-shells. 
Averages are taken over several statistically independent 
realizations of the velocity field. 

Once a stationary state has been reached 10 Million 
tracers are seeded into the flow and integrated according 
to 

X(x ,t) = u(X(x ,t),t), (8) 



where u(X,t) is the velocity obtained from ([6]). 

In order to obtain the velocity at the particle position 
from the grid values we use a tri-cubic interpolation 24 . 
All relevant quantities such as the gradient of velocity are 
stored in intervals of l/7th of the dissipation timescale. 
The main parameters of all simulations are given in Ta- 
ble m 

III. EULERIAN CONDITIONAL STATISTICS 

In this section we present the result on conditioned 
Eulerian increments statistics. We are going to start 
with longitudinal increments, followed by transverse in- 
crements in the subsequent section. 

A. Longitudinal increments 

Following the experimental results obtained by Gagne 
et ali^ - as well as Naert et alJ^, we split our simulation 
domain in subsets f2 e , of fixed rate of energy dissipa- 
tion ei on a line I defined by ((H). On these subsets we 

consider longitudinal velocity increments 5fu in order to 
obtain conditional PDFs P(S^u\ei). The standard (un- 
conditioned) PDFs can be recovered by integrating this 
PDF over all e;. In agreement with the experimental re- 
sults we find for a separation / within the inertial range of 
scales nearly Gaussian statistics for different dissipation 
rates (see Fig. [TJ. Also for different scales I we recover 




FIG. 1. Conditioned PDFs P{5\ u\e{) for different space- 
averaged dissipation rates e; for I = 9377 in comparison to 
a Gaussian distribution, ej = 1 corresponds to the most prob- 
able energy dissipation rate, the others of multiples of this 
rate. All PDFs are normalized to unit variance. 

Gaussianity (see Fig. [2]). The unconditioned PDFs have 
clearly flatter tails. 

As a measure of the deviation from Gaussianity we 
present in Fig. |3] the flatness ((S\ l u) 4 ) / ((S\ l u) 2 ) . The fig- 
ure includes the logarithmic derivative of the third-order 
structure function S^(l) in order to illustrate the inertial 
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TABLE I. Parameters of the numerical simulations. 5Ra = \/15VL/v: Taylor-Reynolds number, it rms : root-mean-square 
velocity, mean kinetic energy dissipation rate, v. kinematic viscosity, dx: grid-spacing, r\ = (v z / e^) 1 ^ 4 : Kolmogorov 
dissipation length scale, t v = (u/e^) 1 ^ 2 : Kolmogorov time scale, L = (2/3£ , ) 3//2 /e^: integral scale, 7l = L/u lms : large-eddy 
turnover time, N 3 : number of collocation points, N p : number of tracer particles. 




FIG. 2. Conditioned PDFs P(5 [ u\ei) for different separations 
I and the most probable value of e, normalized to unit variance 



range via its plateau. The PDFs of the conditioned in- 
crements (Fig. [2]) reach a flatness of approximately three 
throughout the inertial range of scales. 
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FIG. 3. Flatness factors for conditioned velocity increment 
PDFs P(5[ l u\ei), P(6\ l u\Sli), P(S\ l u\Ai) and the uncondi- 
tioned PDF P(5]'u), including the logarithmic derivative of 

the third order structure function Sg ■ The horizontal line 
indicates the flatness of a Gaussian distribution. 

For comparison we conditioned the velocity increments 
on other quantities composed of velocity-gradient tensor 
elements, namely the vorticity ijj — V x u and the longi- 
tudinal gradient I ■ Vw. As for the energy dissipation rate 



we consider spatial averages of the square of vorticity 

1 r l 

9,i = - I dsu\u(x + sl)\ 2 . (9) 
' Jo 

and the square of the longitudinal gradient 

A{' = - / dsv\l- Vu{x + si)\ 2 . (10) 
' Jo 

From Fig. [3] one recognizes that the flatness of 
P(d\ <l u\ei) is closer to a Gaussian distribution over all 

scales than P(S, u\fti) or P(8t u\Ai). It is interesting to 
remark that the integral of the longitudinal gradient over 
I is the longitudinal increment. That the energy dissipa- 
tion rate e; nevertheless works better than this longitudi- 
nal gradient implies that correlations of the form djUidiUj 
with i 7^ j are essential in the condition ([3]) and in the 
RSH ©. 
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FIG. 4. Conditioned PDFs P^j'itjf^) for different separa- 
tions I in comparison to a Gaussian distribution, normalized 
to unit variance 

The flatness of the PDFs conditioned to the scale- 
averaged energy dissipation rate come closest to the 
Gaussian value. However, from Fig. 2] one observes 
that the PDFs conditioned on are also nearly scale- 
invariant but not exactly Gaussian. 

Whether longitudinal increments are conditioned on 
the energy-dissipation or vorticity yields quasi-identical 
results. With respect to the slightly smaller flatness of 
the former one can conclude that longitudinal increment 
statistics is coupled more closely to the scale averaged 
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energy dissipation rate than to the vorticity. This is 
important for diverse models such as the She-Leveque 
model^, where physical reasoning is based on the one 
hand on the energy dissipation rate and the RSH and 
on the other hand on the dimensionality of the coherent 
structures of vorticity. This question is also closely re- 
lated to the issue of different scaling laws for longitudinal 
and transverse structure functions as we will explain in 
the next section. 

We conclude this section on longitudinal increments by 
an examination of the corresponding conditioned struc- 
ture functions S Pj€l . From the scale-invariant PDFs in 
Fig.[3]we expect them to follow the linear K41-scaling law 
p/3 within the inertial range. Indeed, as shown in Fig.[5j 
the conditioned structure functions follow Kolmogorov's 
prediction while the unconditioned higher-order func- 
tions exhibit lower plateaus expressed by a non-zero /i p 
in (P). 



6 i— r 
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FIG. 5. Logarithmic derivative of the conditioned and un- 
conditioned longitudinal velocity structure function of order 
p = 3 and p = 6, lines indicate the K41-prediction 



B. Transverse increments 

In analogy to the RSH for the longitudinal veloc- 
ity increments Chen et al»2& proposed a refined self- 
similarity hypothesis for the transverse velocity incre- 
ments (RSHT). This relation of the scale-averaged square 
of vorticity and the transverse velocity increments reads 

5fu = p 2 (n l l)i, (11) 

where /?2 is a statistical variable independent of I and f2; , 
given by ©. 

Following Chen et al. it is reasonable to look in the 
transverse case at the statistics of velocity increments 
conditioned to f2/, namely the PDFs P(S^u\Ui). They 
are scale-invariant and only slightly flatter than Gaussian 
PDFs (see Fig. [6]). As in the case of longitudinal incre- 
ments one can ask how other conditions such as the en- 
ergy dissipation rate perform compared to e;. The PDFs 
conditioned on e; shown in Fig. [7] are indistinguishable 



from the PDFs conditioned on the vorticity fli. From 
this point of view it is impossible to conclude whether e; 
or fli is the better condition for the transverse fluctua- 
tions. In order to make a more precise statement on the 
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FIG. 6. Conditioned PDFs P(8j L u\€i) for different separations 
I in comparison to a Gaussian distribution, normalized to unit 
variance 

difference between these two conditions it is helpful to 
look at the flatness in Fig. |SJ From this we conclude that 
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FIG. 7. Conditioned PDFs P(S^u\Qi) for different separa- 
tions I in comparison to a Gaussian distribution, normalized 
to unit variance 

conditioning to e; or yields quasi-identical results and 
surprisingly P(Sj L u\ei) are even slightly more Gaussian 
than P(S^u\Qi). The averaged transverse gradient 

1 r l 

= - / dsv\l xWu(x + sl)\ 2 . (12) 
' Jo 

reduces the flatness less than the two other conditions. 
The transverse structure functions are shown in Fig. |H1 
As expected, we find that the conditioned ones follow 
the K41 prediction with the inertial range of scale while 
the high-order unconditioned functions have significantly 
lower plateaus. 
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FIG. 8. Flatness factors of the conditioned PDFs P(5 ; x it|f2;), 
P(Si~u\ei), P(S] L u\A^) and unconditioned P(Sj L u), including 
the third order transverse structure function S3 
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FIG. 9. Logarithmic derivative of the conditioned and uncon- 
ditioned transverse velocity structure function of order p — 3 
and p = 6, lines indicate the Kolmogorov prediction 

IV. LAGRANGIAN CONDITIONAL STATISTICS 

After having computed statistics in the Eulerian frame- 
work we now consider velocity increments ([5]) in the La- 
grangian frame of reference. The Lagrangian analog to 
the RSH might be labeled the Lagrangian refined self- 
similarity hypothesis (LRSH) and reads 

S T v t = /3 L (Te T )^ 2 , (13) 

where the local energy dissipation rate © is averaged 
along a particle trajectory according to 

e T = - [ e(X(x ,t))dt. (14) 
T Jo 

However, one can question whether e T is the correct quan- 
tity appearing in (fT5|) . Benzi et al^i examined this 
relation by means of the assumption of extended self- 
similarity and found that e T rather than the averaged 
square of vorticity 

n T = - f \u(X{x ,t))\dt. 



is the correct quantity in (fl~3|) . Yu et ali^& conditioned 
the velocity increments on a spatially averaged energy 
dissipation rate at one foot-point of the increments. 

In this work we stick to trajectory-averaged condi- 
tions and propose yet another one for Lagrangian in- 
crement statistics. In order to motivate this on di- 
mensional grounds we recall that Eulerian increments 
(ui(l ej) — Ui(0))/l tend to spatial derivatives djUi of the 
velocity field in the limit I — > 0. Those derivatives ap- 
pear in the local energy dissipation rate (j3|). Instead, 
Lagrangian increments (ui(r) —Ui(Q)) / t tend to the fluid- 
particle acceleration in the limit r — > which involve a 
term UjdjUi. We therefore propose to replace (fl~4"l) by 

e^ljdt^d^ + u^r (15) 

hi 

in the LRSH (pjl. 

The calculation of e r , f2 r , and for a given time lag 
r is done by averaging the local quantities over all stored 
points along the particle trajectory. We achieved con- 
verged statistics by taking the average over 10 Million 
particles and several large-eddy turn-over times. 




0.1 1 10 100 

FIG. 10. Flatness factors of the conditioned velocity incre- 
ment PDFs P{8 T Vi\e T ), P(S T Vi\€r) and P(S T Vi\fl T ) as well as 
of the unconditioned PDF P(S T Vi) together with the logarith- 
mic derivative of S% 

In Fig. [10] we compare the flatness of velocity increment 
PDFs conditioned on e T , f2 T , and . We added the loga- 
rithmic derivative of the second order Lagrangian struc- 
ture function S%(1) — ((8 T Vi) 2 ) in order to clarify three 
different ranges of scales: The dissipative scales up to 
r»l, the inertial ones 1 < r < 60, followed by the large 
scales. If we restrict our attention to the inertial range we 
observe that the flatness is most efficiently reduced by e^. 
Also the trajectory integrated energy dissipation rate e T 
diminishes significantly the flatness while the integrated 
vorticity fi T has a negligible effect. This indicates that 

might be a more appropriate condition than e T . 

The corresponding conditioned increment PDFs are la- 
beled P(S T Vi\e T ) and show in Fig. [TTJ The PDF corre- 
sponding to the shortest time-lag considered is reason- 
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FIG. 11. Conditioned PDFs P(8 T Vi\e T ) for different time lags 
t in comparison to a Gaussian distribution and to the K41- 
prediction for the PDF of acceleration, normalized to unit 
variance 

ably well described by the K41-acceleration PDF 16 

P(a) = (a/b)- 5/9 exp[-0.5 (a/6) 8/9 ]/c (16) 

normalized to unit-variance with a — 0.48 and b = 2.72. 
This PDF is the Lagrangian analogon to a Gaussian dis- 
tribution for Eulerian velocity gradients. 

It is important to note that contrarily to the re- 
sults in Eulerian setup the conditioned Lagrangian PDFs 
P(S T Vi\e^) (see Fig. [11]) are still scale-dependent. One 
notes a transition from stretched tails (K41-prediction) 
for short time-lags to Gaussian PDFs (uncorrelated 
statistics) for time lags of the order of the integral time 
scale. This implies that Lagrangian increment statistics 
is 'naturally' scale dependent. 

As can be see from the unconditioned structure func- 
tion in Fig. 1 1 01 Lagrangian structure functions do not 
show a clear scaling law at today accessible Reynolds 
numbers. We therefore refer to relative structure func- 
tions S P (S2)- In the computation of the conditioned 
structure functions we fixed one for all increments 
r. In Fig. Q2] their logarithmic derivatives are shown 
which clearly change under the condition e^. There are 
two major differences between the conditioned and un- 
conditioned functions. The first concerns intermittency: 
The conditioned functions have larger values than the 
unconditioned ones. We observe a value of approxi- 
mately. 1.43 which is close to the K41 prediction of 1.5. 
This implies that intermittency is significantly reduced 
on subsets Q e L . A second feature of Lagrangian incre- 
ment statistics is the so called bottleneck around a few 
r n . It has been attributed to the characteristic trajecto- 
ries (spirals) of tracers in the vicinity of coherent vortex 
filaments. This bottleneck in the local slop is absent once 
velocity increments are conditioned (see again Fig. [T3J) , 
which means that their scaling range is enlarged. Its ori- 
gin is supposed to be in the coexistence of two different 
power-laws. The first related to dissipative effects and 
the second to inertial range physics^,. An insufficient 



separation of dissipative and inertial scales might lead to 
the observed dip in the local slope of structure functions. 
Interestingly, this bottleneck is negligible in the case of 
conditioned structure functions. This implies that it is 
due to a mixture of statistics from different subset ft f L. 




r/r, 



FIG. 12. Logarithmic derivatives of relative Lagrangian ve- 
locity structure function. e T = 1 corresponds to the one with 
the most statistics 



V. CONCLUSION 

This work investigates the statistics of Eulerian and 
Lagrangian velocity increments when conditioned to dif- 
ferent scale-averaged quantities such as the energy dissi- 
pation rate, the square of vorticity or the velocity gradi- 
ent. In the case of Lagrangian increments we propose a 
novel condition dimensionally related to the acceleration 
of fluid elements. 

Considering Eulerian statistics we find that longitudi- 
nal as well as transverse increment PDFs are Gaussian 
shaped with flatness factors close to three when condi- 
tioned to the scale-averaged energy dissipation rate. The 
averaged vorticity produces slightly flatter tails while the 
longitudinal and transverse velocity gradient perform sig- 
nificantly worse. Therefore, there is no preferential link 
of transverse increments and vorticity as of longitudinal 
increments and energy dissipation rate which is impor- 
tant for models of intermittency. Conditional structure 
functions show clear K41-scaling within the inertial range 
of scales. 

Considering Lagrangian statistics we investigated ve- 
locity increments conditioned to trajectory-averaged 
quantities such as the energy dissipation rate, the vor- 
ticity and a novel condition. The latter is motivated by 
dimensional arguments. Conditioning to the dissipation 
rate and to the novel condition yields flatnesses of the in- 
crement PDFs much smaller than without conditioning. 
More precisely, the conditioned PDF of the shortest in- 
crement considered agrees reasonably well with the K41- 
prediction for the PDF of acceleration. Within the iner- 
tial range of scales the flatnesses of PDFs under the novel 
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condition are even smaller than the flatness of PDFs con- 
ditioned on the averaged energy dissipation. 

Conditioned and unconditioned Lagrangian structure 
functions differ significantly. First, conditioning yields 
quasi-K41 scaling exponents. Secondly, the character- 
istics bottleneck of the unconditioned functions at the 
onset of the inertial range disappears once conditioned. 
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